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Abstract 

In this paper we present a version of the general polynomial involutive algorithm 
for computing Janet bases specialized to toric ideals. The relevant data structures 
are Janet trees which provide a very fast search for a Janet divisor. We broach also 
efficiency issues in view of application of the algorithm presented to computation of 
toric ideals. 

1 Introduction 

We consider the problem of computing a Janet basis of a toric ideal in K[x] = . . . , x n ] 
generated by binomials of the form pQ 

T A = { x u - x v | u, v G W l , tt(u) = tt(v), gcd(x u , x v ) = 1 } . 

Here u, v G N n and tt is the semigroup homomorphism 

7r : N n — > Z d , u = {wi, . . . ,u n } — > uia.i H h«„a n 

where a^ G Z d (1 < i < n). 

Given a set of binomials generating a toric ideals, the problem of constructing its Grobner 
basis is usually (except small problems) rather expensive from the computational point of 
view |2j. In practice, for this particular problem, one typically deals with a large number n 
of variables and their degrees. If d is the maximal degree of the initial binomials, then the 
degree of a reduced Grobner basis is bounded by [2j 

d 2 ^ 2 "" 
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But for all that the reduced Grobner basis is also binomial since the binomial structure is 
preserved during the Buchberger algorithm j3J Similarly, the involutive algorithms jB] 
based on the sequential multiplicative reductions of nonmultiplicative prolongations of the 
intermediate polynomials preserve the binomial structure. The output involutive basis which 
is also a Grobner basis though generally redundant. 

Thus, unlike construction of reduced Grobner bases or involutive bases for general polyno- 
mial ideals, the integer arithmetics which may take most of computing time is not important 
for binomial ideals. In this case a fast search of monomial divisors for performing reductions 
of ^-polynomials may become crucial in acceleration of computations. 

Recently [3 |Hj we designed and implemented involutive algorithms specialized to con- 
structing Janet bases of monomial and polynomial ideals. Janet division as well as any 
other involutive division |B] provides uniqueness of an involutive divisor in a polynomial set 
with co-prime leading monomials. This allows one to organize a very fast search for a Janet 
divisor using special data structures for intermediate polynomial sets called Janet trees. 

The main goal of this paper is to discuss the issue of practical efficiency in computing 
Janet bases of toric ideals based on the use of Janet trees. Since one of the most important 
applications of toric ideals is integer programming we shortly describe this application [§| in 
the next section. 

2 Toric Ideals and Integer Programming 

Let A be a matrix of dimension m x n with integer entries and b G Z m , c 6 Z" be vectors. 
The following optimization problem 

min{ c T x | x G N n , Ax = b } 

is called a problem of integer programming. 

We shall assume that c G N n . If there exists vector x satisfying Ax = b, x G N n , then 
the problem of finding minimum of function c T x can be reduced to all kinds of transformation 
of the initial vector state xo using ker(^4). 

The problem of determining ker(A) can be formulated in terms of toric ideals. Indeed, 
every vector u G ker(^4) may be uniquely represented as u = u + — u where both u + and 
u are nonnegative and have disjoint support. Associate symbol Vi with the i-th column of 
matrix A. Then the ideal 

I A = { v" + - v"~ | u+ - u = u G ker(^) } 

associated with ker(^4) is toric. Given the initial solution xo, the optimal solution can be 
found as follows 0: 

1. Construct a basis of the toric ideal X4. 

2. Construct a reduced Grobner basis or an involutive basis of X4 with respect to the 
admissible monomial ordering ^ c generated by vector c. 

3. Reduce monomial v x ° modulo the constructed basis to obtain the optimal solution. 

Therefore, the reduced Grobner basis or any involutive basis of the associated toric ideal I a 
provide an algorithmic tool for solving the problem of integer programming. 
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3 Janet Bases of Toric Ideals 



3.1 Definition of Janet Basis 

In our papers jB] the Grobner bases of special type, called involutive and based on the concept 
of involutive division were introduced. Given a set of coprime monomials and an involutive 
division, any monomial may have at most one involutive divisor in the set. This property of 
the involutive division allows one to design an efficient search for the involutive divisor using 
the method of separative monomials ^U] for a general involutive division or Janet trees [7] 
for Janet division. 

Because of a larger number of variables and unimportance of integer arithmetical oper- 
ations over coefficients of the binomials, the practical complexity of an algorithm for con- 
struction of Grobner or Janet bases is caused by an enormous number of binomials arising 
in computation of the basis. A faster search for divisors may accelerates the computation 
substantially. 

By definition of Janet division [B] (which formalizes the pioneering ideas of Janet [TT] ) 
induced by the order 

x\ y x 2 y . . . y x n (l) 

on x, a polynomial set F is partitioned into the groups labeled by non-negative integers 

[d 1 ,...,d i ] = {f E F | d j = de gj .(lm(/)), l<j<i} 

where degj(-u) denotes the degree of Xi in monomial u and lm(/) denotes the leading mono- 
mial of /. A variable Xj is called (Janet) multiplicative for / G F if i — 1 and 

deg 1 (lm(/)) = max{ de gl (lm(#)) | g E F }, 

or if i > 1, / G [d\, . . . , di-i) and 

degi(lm(/)) = max{ deg f (lm(flf)) | g E [d u . . . , dj_ x ] } . 

If a variable is not multiplicative for / G F, it is called (Janet) nonmultiplicative for /. In 
the latter case we shall write xi E NMj(f,F). u E lm(F) is a Janet divisor of w E M, if 
u | w and monomial w/u contains only multiplicative variables for u. In this case we write 
u \j w. 

Let lm(F) = { lm(/) | / G F }. Then a polynomial set F is called Janet autoreduced if 
each term in every f E F has no Janet divisors among \m(F) \ lm(/). A polynomial h is 
said to be in the Janet normal form modulo F if every term in h has no Janet divisors in 
lm(F). In that follows NFj(f,F) denotes the Janet normal form / modulo F. 

A Janet autoreduced set F is called a Janet basis if 

(V/ E F) (Vx E NMj(f, F)) [ NFj(f ■ x, F) = ] . (2) 

A Janet basis G is called minimal if for any other Janet basis F of the same ideal the 
inclusion lm(G) C \m(F) holds. If both G and F are monic this inclusion implies G C F. A 
Janet basis is a Grobner one, though generally not reduced. However, similarly to a reduced 
Grobner basis, a monic minimal Janet basis is uniquely defined by an ideal and a monomial 
order. In that follows we deal with minimal Janet bases only and omit the word "minimal". 
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3.2 Janet Trees and Search for Janet Divisor 



Consider now a binary Janet tree [7j whose structure reflects the above partition of elements 
in U into the groups which sorted in the degrees of variables within every group. Before 
description of the general structure of Janet trees we explain it in terms of the concrete 
example [7j 

U = {x 2 y, xz, y 2 , yz, z 2 }, (x y y y z) 

and portray it in the form of Janet tree as shown below. In doing so, the monomials in 
set U are assigned to the leaves of the tree. The monomial with increased by one degree of 
the current variable is assigned to the left child whereas the right child points at the next 
variable with respect to chosen ordering. In contrast to Janet tree presented in paper [7j, the 
below tree takes into account sparseness of monomials that is inherent in toric ideals. The 
related information is given in pairs of integers placed in brackets where the first element 
represents the number of current variable and the second one represents its degree. 



.(1,0) 












/ \ (2,0) 






(1,2)\ \ 

\ 2 \ 

\ x y * 


xz 
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\ yz 



Consider now the structure of Janet tree of the general form as a set JT := U{^} of internal 
nodes and leaves which corresponds to a nonempty binomial set. To every element v of 
the tree we shall assign the set of five elements v = {v, d, nd, nv, nb} with the following 
structure: 



var(z/) 


= V 


is 


the index of the current variable 


dg(u) 


= d 


is 


the degree of the current variable 


ndg(i/) 


= nd 


is 


the pointer to the next node in degree 


nvr(z/) 


= nv 


is 


the pointer to the next node in variable 


bnm(z/) 


= bn 


is 


the pointer to binomial 



In the absence of a child we shall assign the value nil to the corresponding pointer. Wherever 
it does not lead to misunderstanding we shall identify the pointers nd and nv with the nodes 
they point out. To the root of JT we assign u with var(i/ ) = 1 in accordance with labeling 
© and dg(v ) = 0. 
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The internal nodes and leaves of tree JT are characterized by the states: 

Internal node: {{nv ^ nil A v < var(nu)) V (nd 7^ nil A d < dg(nd)) 
A bn = nil 

Leaf: nv = nil And = nil A bn 7^ nil A d = dg(lm(6n)). 

For a fast search for Janet divisor in the given tree one can use the following algorithm 
J-divisor which is an adaptation to the above structure of Janet tree of the algorithm 
described in j7]. 



Algorithm: J-divisor( JT, w) 



Input: JT, a Janet tree; w, monomial 




Output: bn, a binomial such that lm(&n) 






or nil, otherwise 




1 


v := u 




2 


while deg vai(l/) {w) > dg(u) do 




3 


while ndg(i/) and deg var(ndg(!y)) (w;) > 


dg(ndg(z/)) do 


4 


v := ndg(z/) 




5 


od 




6 


if nvr(z/) then 




7 


v := nvr(z/) 




8 


elif ndg(z/) then 




9 


return nil 




10 


else 




11 


return bnm(^) 




12 


fi 




13 


od 




14 


return nil 





Apparently, the next theorem formulated and proved in [7] is valid for the adapted algorithm 
as well. 

Theorem . Let d be the maximal total degree of the leading monomials of binomials in n 
variables which constitute the finite set U . Then the complexity bound of the algorithm 
J — divisor and the binary search algorithm is given by 

ij-divisor = 0(d + n), 

^BinarySearch = 0(n({d + u) log(rf + u) -Tllog(n) - dlog(d))). 

Thus, the complexity bound for the search of Janet divisor is 0(n + d) where n is the number 
of variables and d is the maximal degree of the leading monomials in the binomial basis. 
Since this bound is even lower than that for the binary search algorithm, one can expect 
that the involutive completion of binomial ideals may be faster than the reduced Grobner 
basis completion. 
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3.3 Algorithms for Binomial Janet Bases 

Given the generating binomial set F of a toric ideal Z_4, the following algorithm Binomial- 
JanetBasis which is a special form of the general polynomial algorithm [HJ El constructs a 
Janet basis of X4. 



Algorithm: BinomialJanetBasis(F, -<) 



Input: Fel\{0}, a finite binomial set; -<, an admissible 




ordering 


Output: G, a Janet basis of the ideal generated by F 


1 


choose / G F with the lowest lm(/) w.r.t. y 


2 


T:={/,lm(/),0} 


3 


Q:={ {g,lm(g),0} \qeF\{f} } 


4 


Q :=JanetReduce(Q, T) 


5 


while Q ^ do 


6 


choose p G Q with the lowest lm(bin(p)) w.r.t. >- 


7 


Q-=Q\ M 


8 


if hn(bin(p)) = anc(p) then 


9 


for all { r G T \ Zm(bin(r)) y lm(bin(p)) } do 


10 


Q:=QU{r}; T := T \ {r} 


11 


od 


12 


p:= NFj(bm(p),T) 


13 


fi 


14 


T := T U {p} 


15 


for all q G T and x G A^Mj(bin(g), T) \ nmp(g) do 


16 


Q := Q U { {bin(g) • x, anc(g), 0} } 


17 


nmp(g) := nmp(g) H iVMj (bin(g), T) U {x} 


18 


od 


19 


Q :=JanetReduce(<5, T) 


20 


od 


21 


return G := { bin(/) | / G T } 



As well as in |Hj to apply the involutive criteria and avoid repeated prolongations we shall 
endow with every binomial / G F the triple structure 

p = {f, u, vars} 

such that 

bin(p) = / is binomial itself, 

anc(6) = u is the leading monomial of a binomial ancestor of / in F 
nmp(p) = vars is a (possible empty) subset of variables. 

Here the ancestor of / is a polynomial g G F with u = lm(g) and such that u \ \m(p). 
Moreover, if deg(w) < deg(lm(p)), then every variable occurring in the monomial \m(p)/u 
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is nonmultiplicative for g. Besides, for the ancestor g the equality anc(g) = lm(g) must 
hold. These conditions mean that polynomial p was obtained in the course of the below 
algorithm BinomialJanetBasis from g by a sequence of nonmultiplicative prolongations. 
This tracking of the history in the algorithm allows one to use the involutive analogues of 
Buchberger's criteria to avoid unnecessary reductions. 

The set vars contains those nonmultiplicative variables that have been already used in 
the algorithm for construction of nonmultiplicative prolongations. This set serves to prevent 
repeated prolongations. 

In order to provide minimality of the output Janet basis we separate IE] the whole 
polynomial data into two subsets which are contained in sets T and Q. Set T is a part of 
the intermediate binomial basis. Another part of the intermediate basis is contained in set 
Q together with all the nonmultiplicative prolongations of polynomials in T which must be 
examined in accordance to the above definition (0) of Janet bases. In so doing, after every 
insertion of a new element p in T all elements r G T such that 

lm(bin(r)) >- lm(bin(p)) 

are moved from T to Q as the for-loop 6-11 in algorithm BinomialJanetBasis does. Such 
a displacement guaranties that the output basis is minimal [U]. 

It should also be noted that for any triple p G T the set vars must always be a subset of 
the set of nomultiplicative variables for binfjo) 

vars C NMj(bm(p),T) . (3) 

In the description of algorithm JanetBinomialBases we use the contractions: 

NMj(bm(p),T) = NMj{bm(p), {bin(/) | / G T}) , 
NFj{bm(p),T) ee NFj{bm(p), {bin(/) | / G T}) , 

The insertion of a new polynomial in T may generate new nonmultiplicative prolongations 
of elements in T which are added to Q in line 16. To avoid repeated prolongations the set 
nmp(g) of Janet nonmultiplicative variables for g has been used to construct its prolongations 
is enlarged with x in line 17. The intersection placed in this line preserves the condition (jSJ). 

The subalgorithms JanetReduce and NFj perform Janet reduction of polynomials in 
Q modulo polynomials in T and presented below. In addition to reductions in lines 4 and 
19, the Janet normal form computation is placed in line 12. This is because the replacement 
of elements from T to Q may lead to the tail reducibility of the binomial in p. Such a 
reducibility may be caused by converting of some nonmultiplicative variables for binomials 
in T into multiplicative due to the replacement. 

In subalgorithm JanetReduce computation of the Janet normal form h is done in line 6 
for every binomial bin(p) in T. If h is nonzero, then line 8 checks if lm(bin(j»)) was subjected 
by reduction. If the reduction took place lm(/i) cannot be multiple of any monomial in the 
set {lm(bin(g)) | g G T} [6 . Therefore, one has to insert the triple with h in the output set 
Q as shown in line 9 as h cannot have ancestors among polynomials in T and one must also 
examine all nonmultiplicative prolongations of h. If lm(bin(p)) is Janet irreducible modulo 
{lm(bin(g)) | g G T }, then the triple {h, anc(p), nmp(p)} is added to Q in line 11. 
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Algorithm: JanetReduce(Q, T) 



Input: Q and T, sets of triples 


Output: Q whose polynomials are Janet head reduced modulo T 


1 


S := Q 


2 


Q ■= 


3 


while S ^ do 


4 


choose p E S 


5 


S:=S\{p} 


6 


h := NFj(jo,T) 


7 


it ai 7^ U then 


8 


if lm(bin(p)) ^ \m(h) then 


9 


g:=gu{/i,lm(/i),0} 


10 


else 


11 


Q := QU {h, anc(p), nmp(p)} 


12 


fi 


13 


fi 


14 


od 


15 


return Q 



Subalgorithm NFj(p, T) performs the Janet reduction of a binomial g = bin(p) modulo 
polynomial set in T: 

Algorithm: NFj(/,T) 



Input: / = {bin(/), anc(/), nmp(/)}, a triple; T, a set of triples 


Output: h = NFj(bin(f),T), the Janet normal form of the 




binomial in / modulo binomial set in T 


1 


G := {bin(s) | g E T} 


2 


h := bin(/) 


3 


if lm(h) is Janet reducible modulo G then 


4 


choose g G T such that lm(bin(g)) \j lm(h) 


5 


if \m(h) 7^ anc(/) and 




CriterionI(/, g) or CriterionII(/, g) then 


6 


return 


7 


fi 


8 


else 


9 


while h 7^ and /i has a term t Janet reducible modulo G do 


10 


choose q E G such that lm(g) j t 


11 


h := h — q ■ t / lm(g) 


12 


od 


13 


fi 


14 


return h 



For the head reducible input binomial bin(/) the two criteria are verified in line 5: 
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• CriterionI(/, g) is true iff anc(/) • anc(g) | lm(bin(/)). 

• CriterionII(/, g) is true iff deg(lcm(anc(/) • anc(g))) < deg(lm(bin(/)). 

These criteria are the Buchberger criteria [12] adapted to the involutive completion proce- 
dure. If any of the two criteria is true, then NF(bin(f), T) = (Hi- 
lt should be noted that the Janet normal form is uniquely defined and, hence, uniquely 
computed by the above subalgorithm. This uniqueness hold because of the uniqueness of 
a Janet divisor among the leading terms of binomials in T at every step of intermediate 
computations [H]. 



4 Examples 

As we emphasized in Sect. 3.1, in the course of involutive completion of the initial binomial 
generators for a toric ideal the reduction can be performed very fast due to the fast search 
for a Janet divisor, This fast search is provided by the use of the Janet tree structures for 
intermediate binomial set. Our computer experiments with C/C++ codes implementing 
polynomial algorithms for Janet bases jH] perfectly strengthen this theoretical fact. In par- 
ticular this fast reduction in addition to suppressing swell of intermediate integer coefficients 
results in high computational speed observed for the benchmark collection used for testing 
Grobner bases software jSj. These benchmarks, however, are not very "sparse" with respect 
to degrees of variables occurring in the generating set. By contrast, the generating binomial 
sets for toric ideals especially for those arising in integer programming problem are usually 
highly sparse. This may lead to much larger cardinality of a Janet basis than that of the 
reduced Grobner basis and thereby annihilate the advantages of involutive reduction. 
Consider the example taken from ^3] 

* / . f HP ry ry „ ry ry 1 ry^^ ?y ff 1 4 „ 2 /-y» 39 < -y l 25~,14 "1 

-L' J\ — ^ j q J- \ Ju 2 <*j 3 4 -L 5 3 1 ^4 j 1 2 3 I ' 

Our C++ package |H] generates the degree-reverse-lexicographical Janet basis of X4 with 
7769 binomials whose sorting with respect to the ordering chosen gives 

f S 281 280 fil 2 221 27Q 2 281 280 1 1 

J ry ry ^ rp ry ry ry ry ry^ ry ry ry ry ' " ry ry ry ry ^ ry ry ry ry ry ry I L 

%k> ^ ]^ 3 4 jw^ju^ 5 tX- 2 3 4 ju\<ju<£ - iX.' q iX. ^ tX.' 3 *^ 4 2 5***5 •x-' Q «•<-•' 1 2 ^ 3 4 j 

where we explicitly show only three highest ranking binomials and the lowest one. The 
computing time on a Pentium III 700 Mhz based PC running under RedHat Linux 6.2 is 6 
seconds that is noticeably larger than the running time for direct computation of the reduced 
Grobner basis which contains 19 binomials only: 



r 2 281 280 281 280 2 221 218 2 219 220 

\ XqX\ X3X4 — X2 , X2 — X1X4 , X0X3 X4 — X\X2 , X\ X2 — X3X4 , 

~ n, 3„ 161 4 156 5_ 157 2 160 4_ 101 7_ 94 8_ 95 3_ 100 

X0X3 X4 — Xi X2 , xi X2 — X3 X4 , X0X3 X4 — xi X2 , xi X2 — X3 X4 , 

X Xi 4 X 4 61 - X 2 61 , X 2 62 X 3 - Xi 3 X 4 6 °, X X 3 5 X 4 41 - X! 10 X2 32 , Xi U X 2 33 - X 3 4 X 4 40 , 

~ ~ 26„ 15 _ _ 38 39 25„ 14 15_ 21 28_ 4 29 5 14 20 

X0X2 X3 X4 — Xi , Xi — X2 X3 , X0X1 X4 — X2 X3 , X2 X3 — X\ X4 , 

10 21 24 3 25 4 9 20 1 1 

X0X3 X 4 — Xi X2 , Xi X2 — X3 X4 , X0X1X2X3X4 — 1 j • 

Accordingly, such a computer algebra system as Singular [T7| needs much less than 1 second 
to compute this Grobner basis on the same computer. 



9 



Having ascertained this drawback of the involutive method with respect to the Grobner 
basis one in computing toric ideals we designed another algorithmic approach to computing 
Grobner bases ^3]. This approach preserves the Janet-like tree structure and uniqueness of 
a divisor though underlying division is not involutive since it does not satisfy the axioms 
in jH] . On the other hand the resulting bases unlike Janet bases are often reduced as Grobner 
bases and their cardinality is always less or equal to the cardinality of Janet bases. For toric 
ideals the new bases are much more compact then Janet bases. We have not implemented 
yet the new algorithm and so we demonstrate the compactness of its output in comparison 
with algorithm Binomial JanetBasis by the following simple example taken from [2]: 

Aa = { x — y z,x w — y ,x y — zw } . 

The reduced Grobner basis and Janet basis of this toric ideal for the degree-reverse-lexico- 
graphic order induced by x>~y>~z>~w are 

r 7 2 4 4 2 1 

{ x — y z,x w — y ,x y — zw, y — xzw } 

and 

r 7 2 fi 3 fi 2 3 5 2 24 3 25 3 

{ x — y z,x y — x zw, x w — x y ,x y — x zw, x y — x zw ,x w — xy , 
x 4 y — xzw, x 2 zw 2 — xy 4 , x 4 w — y 3 , x 3 y — zw, y 4 — xzw 2 } , 

respectively. Their cardinalities are 4 and 11. The new basis contains 5 elements 

{ x 7 — y 2 z, x 4 y — xzw, x 4 w — y 3 , x 3 y — zw, y 4 — xzw 2 } 

and contains only single extra element in comparison with the reduced Grobner basis. 

It should be noted that there are also a number of other efficient algorithms computing 
Grobner bases of toric ideals (see, for example, which are differ greatly from just 

completion of a generating binomial set to a Grobner basis. After implementation of our 
new algorithm we are planning to run the underlying code for collection of large examples 
given in Pfl|2] and other references. 
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